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Hierarchical formation of bulgeless galaxies: 
Why outflows have low angular momentum. 
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ABSTRACT 

Using high resolution, fully cosmological smoothed particle hydro-dynamical sim- 
ulations of dwarf galaxies in a A cold dark matter ( ACDM) Universe, we show how 
high redshift gas outflows can modify the baryon angular momentum distribution 
and allow pure disc galaxies to form. We outline how galactic outflows preferentially 
remove low angular momentum baryons, due a combination of a) star formation peak- 
ing at high redshift in shallow dark matter potentials, an epoch when accreted gas 
has relatively low angular momentum, b) the existence of an extended reservoir of 
high angular momentum gas in the outer disc to provide material for prolonged SF 
at later times and c) the tendency for outflows to follow the path of least resistance 
which is perpendicular to the disc. We also show that outflows are enhanced during 
mergers, thus expelling much of the gas which has lost its angular momentum during 
these events, and preventing the formation of "classical" , merger driven bulges in low 
mass, field galaxies. Stars formed prior to such mergers form a low surface brightness, 
extended stellar halo component, similar to those detected in nearby dwarfs. 

Key words: galaxy evolution-galaxy formation-cosmology. 



1 INTRODUCTION 

The angular momentum of disc galaxies is thought 
to originate from tidal torques imparted by surround- 
ing structures in the exp anding Universe l|Pee bles 1969; 
iBarnes fc Efstathioulll987l ). prior to proto-galactic collapse. 
Disc galaxies will form if gas gains a similar amount of an- 
gular momentum as the dark matter in this process, and if 
this angular momentum is subsequently retained as the gas 
cools to the centres of dark m atter halos and settles into a 
disc ijFall fc Efstathioulll980h . Gas subsequently fragments 
and forms stars. Within cold dark matter (CDM) cosmol- 
ogy this picture of the origin of galactic angular momentum 
is believed to hold even as galaxies are built hierarchically, 
with the added assumption that mergers result in building 
bulge components as discs are d estroyed by violent relax- 
ation l|Barnes fc Hernauistlll99"^ ). 

However, detailed modelling has highlighted problems 
with this scenario. In simulations, gas cools to the centre 
of the earliest collapsing structures making them very cen- 



trally concentrated. Dynamical friction occurring during the 
mergers of such systems results in significant loss of angu- 
lar m omentum (jNavarro fc Steinmetzfeood : iMaller fc Dekell 
2002). This problem has largely been overcome in simula- 
tions by i ncreasingly effective recipes for feedback from su- 
pernovae jThacker fc Couchmanl |2000| ; IsTinson et al.ll2006l ; 
iGovernato et al.l 120071 ) . and increased resolut ion that de- 
crease spurious angul ar momentum losses (jMaver et al.l 
|200S| ; ISales et alj|201(ih . Despite this encouraging progress, 
decreasing the impact of dynamic friction does not solve 
all the problems of galaxy formation which relate to angu- 
lar momentum within hierarchical structure formation. In 
addition, the angular momentum distribution of baryons 
in real galaxies differs signifi cantly from the distribu- 
tion within dark matter halos l|van den Bosch et al.l 1200 ll ; 
Dut ton fc van den Boschl l2009). Even if angular momentum 
were fully co nserved, CDM halo s have a low angular mo- 
mentum tail |Bullock et al.ll2~00l[ ). Recent simulations show 
that galaxies formed in a CDM Universe have realistic disc 
sizes I Brooks et al.|[201ol ) , however the close coupling of dark 



2 C. B. Brook et al. 



matter and gas while acquiring their angular momentum 
implies that baryons should also have a substantial low an- 
gular momentum component, which is generally believed to 
be reflected in the bulge components l|van den Boschll200ll ; 
IStinson et alj201oh . Yet many real galaxies have no bulge; ~ 
70% of galaxies with stellar mass o f < 10 9 Mq have a Sersic 
index of less than 1.5 (|Duttonl2009l ) and many more massive 
galaxies do not have classical bulges (|Kormendv et al]|2010l ; 
iFisher fc Drorvl feoiO). Pure disc structures reflect an under- 
lying angular momentum distribution which is seemingly at 
odds with theorie s of galaxy formation within the accepte d 
ACDM paradigm (jBullock et al.ll200ll :l van den Bosch 2001). 
Furthermore, in a hierarchical CDM Universe, mergers drive 
gas toward the centre of DM halos, making the formation of 
low angular momentum bulge stars an apparently inevitable 
consequence. The persistence and degree of these problems 
means that the manner in which disc galaxies attain their 
angular momentum distribution has been considered one of 
the major long standing challenges to the ACDM paradigm. 

The removal of low angular momentum gas in galactic 
outflows is one viable solution to explain the discrepancy 
in angular m omentum distributi o ns between dark m atter 
and baryons jPekel fc Silk! 1 19861 ; iBinnev et all 1200 ll ) and 
the main focus of this work. Several lines of evidence sug- 
gest that supernova driven galactic winds are able to ex- 
pel large amounts of gas from galaxies dur i ng the galaxy 
form ation process (|Mathews fc Bakerl [l97ll ; IVeilleux et al.l 
120051 ). Galactic winds have been observe d, and extensively 
studi e d in local star forming galaxies (e.g. | Lvnds fc Sandagd 
19631; lAxon fc Taylorl Il97si; Martini Il999t lOhvamal 12003 : 
Martin et al. 2005; Rupkc et al. 2005r i. or inferred from the 



metallicit y distribution of stellar populations in low mass 
galaxies (|Kirbv et all l201ll ). Outflows are expected to be 
more c ommon at high reds hift where star formation is more 
active (|Madau et al.|[l996l ). Indeed, outflows have been de- 
tected at high redshift from Lyman-break selected galaxies 
JPettini et al.ll 19981: ISimcoe et all [20021 : IShaplev et atll2003l : 



lAdelberger et al.| 2005), as well as gravitationally lensed 
galaxies ( Pettini et al.l 2000! ) . 

We recently simulated dwarf galaxies which match a 
wide range of observed dwarf galaxy properties, includ- 
ing the absence of a bulge and cored DM distribution 
that ca uses a linearly rising ro tation curve with no inner 
"peak" (|Governato et all 120101 ; G10 hereafter). The simu- 
lated galaxies lie simultaneously on t he Tully-Fisher rela - 
tion, and the size-luminosity relation (Bro oks et al.| [2010). 
Using techniques which closely mimic observations, the 
galaxies have been shown to have mass distributions and 
stellar masses w hich are consist ent with galaxies in the 
THINGS survey l|Oh et al l (|2010l ). Figure 5) of similar ro- 
tational velocity and fall on the baryo nic Tully-Fisher re- 
lation, as derived by iMcGaughl (|2005l ). Most importantly, 
the simulated dwarfs also match the distribution of angular 
momentum of stars observed in pure disc galaxies, differing 
significantly from that of the dark matter halo in which they 
are embedded (see Figure 4 of G10, as well as Figure [8] of 
this paper). These properties indicate that our simulations 
can provide unique insights into the acquisition and reten- 
tion of angular momentum of disc galaxies. Here we pro- 
vide an overview of the various interlinked processes which 
determine the angular momentum distribution of stars in 
our simulated bulgeless disc galaxy. In particular we outline 



why the significant amount of gas which is expelled from the 
galaxy is primarily low angular momentum material. 

In Section [5] we review the code and the simulation 
initial conditions. We present the properties of our simu- 
lated galaxy in Section \2. 11 highlighting resolution issues as 
well as the effects of allowing gas to cool to dense regions 
before forming stars. We plot star formation rates and 
outflow rates, as well as their ratio (the "mass loading") 
in Section We show that the outflows have a strong 
bias toward low angular momentum in Section [3] We 
outline the reasons for outflows having preferentially low 
angular momentum material: early accreted material has 
low angular momentum, yet needs to escape a relatively low 
potential well (Section 14. 1[) ; the existence of an extended 
reservoir of high angular momentum gas (Section 14. 2p : 
with outflows being perpendicular to the disc (Section 14, 3p . 
We also show that the starburst triggered in the late 
merger expels the bulk of the gas which lost its angular 
momentum during this event (Section I4.4|) . Our summary 
and discussion follows in Section [5] 



2 THE SIMULATION. 

At z = the virial mass of the halo in the dark mat- 
ter only run of our simulated galaxy is 3.7 x 10 10 Mq 
(the virial mass is measured within the virial radius R v i r , 
the radius enclosing an overdensity of 100 times the cos- 
mological critical density). The halo was selected within 
a large scale, low resolution, simulation run in a concor- 
dance, flat, A-dominated cosmology: Qp = 0.24, A=0.76 , 
h = 0.73, erg = .77, and n b = 0.042 (|Spergel et al.ll2003l : 
I Verde et al.ll2003T ). We evolved the simulation using the fully 
parallel, N-body+smoothed particle hydro-dynamics (SPH) 
code GASOLINE which self consistently models gas hydro- 
dynamics and cooling, star formation, stellar mass loss, en- 
ergy feedback from supernovae, metal enrichment and gas 
diffusion together with structure formation within theCDM 
model. A detailed description of the code and of the sim- 
ulations is available in the literature ( Wadslev et all 120041 : 
iGovernato et al.ll2010l : IShen et al.ll201ol ). The version of the 
code used in this study includes radiative cooling and ac- 
counts for the effect of a uniform background radiation field 
on the ionization and excitation state of the gas. The cosmic 
ultraviolet bac kground is implemented using the Haardt- 
Madau model |Haardt fc Madaulll996l ). we use a standard 
cooling function for a primordial mixture of atomic hydro- 
gen and helium at high gas temperatu res, but we include 
the effect of metal coolin g for T < 10 4 (|Bromm et al.ll200ll : 
iMashchenko et al"1 l2008). The possible effects of including 
H2 metal line cooling and SF linked to H2 abundance and 
some preliminary results are briefly discussed in §5. 

Star formation occurs when cold g as reaches a given 
threshold density (|Stinson et all 120061 ). which is depen- 
dent on the vo lume of the star forming regions that 
can be resolved ( Saitoh et al. 200§|; iRobertson fc Kravtsovl 
l200Sl : iTasker fc Brvanll2008l ). We used 100 amu/cm 3 in the 
main run analysed here. The adoption of a higher density 
threshold for star formation (most previous studies used 



O.lamu/cm ) has the effect of limiting SF to gas clumps 
similar in mass to real star forming regions and of increas- 
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ing drastically the amount of energy per unit mass releasi 
into the gas directly affected by supernova feedback. Tl 
higher threshold is justified by the high mass and spati 
resolution of this run that allow individual star forming i 
gions to be resolved. Star formation then proceeds at a ra r- 
proportional to pl'as, i.e. locally enforcing a Schmidt la -c » 
The adopted feedback scheme is implemented by releasi] 
thermal energy fr om supernovae into the gas surroundii -3* 
each star particle ijStinson et al.ll2~006l l. The energy relea 
rate is tied to the time of formation of each particle (whk_200 
effectively ages as a single stellar population with a Krou] 
IMF). To model the effect of feedback at unresolved scak 
the affected gas has its cooling shut off for a time scale pr 
portional to the Sedov solution of the blastwave equatio 
which is set by the local density and temperature of t 
gas and the amount of energy involved. In the main run d 
scribed in this study this translates into regions of ~ 150 
to 400 pc in radius being heated by feedback from multiple, 
overlapping supernovae, and having their cooling shut off 
for typically 5-10 million years. However, even during high z 
starbursts only a few percent of the gas in the disc plane is 
heated by supernova feedback to temperatures > 40,000 K. 
The effect of feedback i s to regulate star fo rmation efficiency 
as a function of mass ijBrooks et al.ll2007t l. Only two other 
parameters are needed, the star formation efficiency (eSF = 
0.1) and the fraction of supernova energy coupled to the in- 
ter stellar medium (ISM) (eSN = 0.4). Si milar values have 
been used in previous works by this group |Governato et all 
l2007h . However, here we slightly increased eSF from 0.05 to 
0.1 to ensure a better normalization to the local Schmidt 
law. We explored values of eSN as high as 0.6 and cooling 
shutoff times changing by a factor of a few to verify that 
results are robust to small changes in the description of SF. 

2.1 The Galaxy: Resolution and Cosmic Variance 

The simulation presented here is part of an ongoing cam- 
paign to study the formation and evolution of galaxies in a 
ACDM cosmology. Recent work has highlighted the role of 
numerical resolution in driving the structural properties of 
simulated gala xies and specifically their internal m ass distri- 
bution (jMaver et al.ll200Sl : IGovernato et alJl2010h . We anal- 
ysed the simulated galaxy used in this study at different res- 
olutions and showed that our key results are retained at high 
resolution regimes, which have spatial resolution less than 
~ 200 pc and gas particle mass less than a few 1O 4 M (G10). 
Yet we emphasise the importance of resolving star forming 
regions by allowing gas to collapse to high density regions 
during these merger events. For example a high SF thresh- 
old DG1MR forms about 3x10 s M in stars. This rela- 
tively low star formation efficiency leads to an M s t ars /Mh a ; 
ratio of ~ 1/200, which is close to what is measured for 
nearb y small galaxies with resolved kinematics and photom- 
etry l|0h et alJfeoiCl ). On the other hand a 0.1 amu/cm 3 
threshold allows eight times more stars to form (see G10, 
online material). Adopting high resolution significantly de- 
creases SF also in halos with mass > 10 11 Mq (Governato et 
al, in prep, Mayer et al in prep ) , reducing the overproduc- 
tion of stars (|Brooks et alj|20ich and bringing them cl ose to 
estimates for more massive galaxies (|Guo et alJl20ld ) that 
predict M 3tars /Mh a ( ~ 0.04-0.05. The importance of resolv- 
ing star forming regions was also demonstrated by the results 
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Figure 1. A type of merger tree. We plot the birth position within 
the simulation box of every star which ends up in the final galaxy 
at z=0. We are tracing the birth of stars as the galaxy moves 
through the simulation box. The axis are the x and y co-ordinates 
of the simulation box, shown in co-moving kpc. The colors are the 
birth time of the stars, where 2 = corresponds to 13.7 Gyrs.. We 
see in this plot the merger events which involve satellites which 
are large enough to host stars. Two significant merger events in 
particular can be seen at ~ 2.5 Gyrs and ~ 7 Gyrs. For these 
events, we note the virial mass of the satellites (Satl and Sat2) 
prior to the merger, as well as the mass of the main progenitor 
at corresponding times. The passages of Sat2 around the central 
galaxy and its extended merging period is also evident. 
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Figure 2. The star formation rate plotted as a function of time. 
We show 200 Myr bins in red, highlighting the bursty nature 
of the star formation. The black line shows 1 Gyr bins, which 
is shown because that matches the time resolution information 
that we have for outflows, which are less well defined than star 
formation. The peak in star formation at ~ 7Gyrs is associated 
with the last major merger. 



from a constrained merger simulation done bv lTevssier et al] 
(2010). The more localised star formation within an inho- 
mogeneous inter stellar medium is crucial in our model in 
terms of allowing large scale outflows to be generated (G10). 

From the simulation we generated frequent outputs 
in order to determine outflow rates and the properties 
of the gas prior to being blown out. To this end, we use 
DG1MR from G10, which has outputs every ~ 500 Myrs, 
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Figure 3. The outflow rate of gas from the galaxy as a func- 
tion of time. The black line shows the outflows from the central 
galaxy. The red dashed line adds outflows from the two most 
significant satellites accreted during the galaxy's evolution. The 
peak outflows occur during the merger events at ~ 2.5 Gyrs and 
~ 7 Gyrs. 
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Figure 4. The evolution of the ratio of the mass of expelled gas 
to the mass of stars formed, showing a 'mass loading' from the 
central galaxy (square symbols) of greater than one for almost the 
entire simulation, and a rise associated with a late major merger 
event. The crosses show the mass loading of the major satellite 
which merges at Z ~ 1, which also shows significant outflows. 

and has an effective resolution 2304 3 particles in a 25 Mpc 
box, with gravitation softening of HOpc and initial stellar 
mass of ~ 17OOM0. The morphology (that of a bulgeless 
dwarf irregular galaxy), star formation rate, final stellar 
and baryon fraction, Mhj/Lb, rotation velocity, magnitude 
and color show little change between DG1 and DG1MR. 
These properties contrast significantly compared to the 
same simulation run with a low threshold for star formation 
(DG1LT), which does not generate large scale outflows. 
Details of the properties of the galaxies DG1, DG1MR and 
DG1LT using techniques which mimic observations were 
presented in G10. DG1MR has M, = -16.9, g-r = 0.54, 
M star = 5.3 x 1O 8 M0, scalelength R s = 0.9kpc, rotational 
velocity V ro t = 55kms _1 , and ratio of neutral hydrogen to 
B-band luminosity Mhi/Lb = 1.0. 

Bulges are commonly asserted to be formed in merg- 
ers, with galaxies which have more mergers, or mergers 
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Figure 5. The cumulative mass of outflows as a function of cir- 
cular velocity. Significant outflows occur prior to the simulated 
galaxy attaining all its material. The top axis indicates the cor- 
responding redshifts. 

with larger mass ratios, assumed to have larger b ulges 
(e.g. I Barnes fc Herncmistl Il996l ; iHopkins et ail |2010|) . In 
CDM cosmologies, merger histories of low mass galaxies are 
equally as diverse and rich as high mass galaxies, yet it is 
only i n low mass fie ld galaxies where bulgeless discs are com- 
mon (jDuttonll2009[) (dwarf satellite s often have spheroidal 
morphologies e.g. iMaver et al.ll20oih . To validate the CDM 
model, simulations should form low mass field galaxies with 
low B/D ratios and rising rotation curves starting from most 
realizations of cosmological initial conditions. It is important 
to highlight that the success of our analysed galaxy is not a 
result of a particular, special merging history. We showed in 
G10 the results of two simulations which had very different 
merging histories, highlighting that we have not selected a 
particularly quiescent simulated galaxy. In Figure [T] we plot 
a type of merger tree for the simulation analysed in this 
study, which shows the birth position within the simulation 
box of every star which ends up in the final galaxy at z=0. 
We are tracing the birth of stars as the galaxy falls through 
the simulation box. The axes are the x and y co-ordinates of 
the simulation box, shown in co-moving kpc. The colors are 
the birth time of stars. We see in this plot the merger events 
which involve satellites which are large enough to host stars. 
Two significant merger events in particular can be seen at 
~ 2.5 Gyrs and ~ 7 Gyrs. Prior to merging, Satl has a virial 
mass of 1.4e9Me, S as mass of 2.Oe8M0 and stellar mass 
of 3.7e6Mo, with the central galaxy having virial mass of 
4.6e9Mo, gas mass of 2.6e8M0 and stellar mass of 4.1e7M0 
at a corresponding time. Also evident are the passages of 
Sat2 around the central galaxy and hence its merging period 
which extends for ~ 1.6 Gyrs between ~ 5.2Grys (z = 1.2) 
and final coalescence at ~ 6.8 Gyrs (z — 0.8). Prior to merg- 
ing, Sat2 has a virial mass of 4.9e9M0 , gas mass of 3.4e8M0 
and stellar mass of 3.8e7M0, while the central galaxy has 
virial mass of 1.4elOM0, gas mass of 1.3e9M0 and stellar 
mass of l.SeSM© at a corresponding time. 



3 STAR FORMATION AND OUTFLOWS 

In Figure [2] we plot the star formation rate as a function 
of time for our simulated galaxy. In red, we show the star 
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Figure 6. The evolution of the ratio of the mass of expelled gas 
to the total virial mass of the galaxy. The ratio peaks at early 
times, with the subsequent steady decline interrupted by a rise 
which occurs during the last major merger event. 



formation using 200 Myr bins. The bursty nature of star 
formation is apparent. The black line shows 1 Gyr bins, to 
allow a better match to the time resolution we have for our 
outflows in subsequent plots. Star formation shows bursts, 
some of which are associated with mergers, with a high star 
formation rate at ~ 2.5 Gyrs associated with Satl, and then 
peaks again at ~ 7 Gyrs. This final increase in the star for- 
mation rate is associated with the last major merger (Sat2), 
which has a mass ratio of ~ 3 : 1. Subsequent to the end of 
this merger, the star formation rate decreases significantly, 
and is quite low for the final ~ 4 Gyrs but continues to be 
bursty. 

The outflow rate, in yr~ , is shown in Figure [3] The 
black line shows the outflows from the central galaxy. The 
red dashed line adds outflows from the two most significant 
satellites accreted during the galaxy's evolution, marked 
Satl and Sat2 in Figure [1] Outflows are calculated in a very 
simple manner: any gas which is bound and within 5R* at 
any output timestep, but is outside R V i r of the final galaxy 
at z = 0, is assumed to have been expelled, where R, is the 
radial distance to the furthest young star (age < 200 Myrs) 
from the centre of the galaxy. Typically, 5it* is ~ 0.5R V i r - 
Our results are not sensitive to this choice, so long as the ra- 
dius is well beyond the star forming region which means that 
we do not miss gas which is being expelled, and is not too 
close to R V ir, which can result in gas which is not bound be- 
ing falsely identified as outflows. We only include gas which 
does not return to the galaxy at later times in our definition 
of outflows in this study. Using 28 outputs which are ~ 500 
Myrs apart, each outflow particle is traced through every 
output to find the maximum "jump" in radial distance from 
the centre. The jump in radius is invariably associated with 
a temperature that shows that the gas particle was heated 
by supernova energy. The time of the radial and temperature 
jump determines the outflow time for the particle. 

We note that the degree and direction of out- 
flows in our simulation is n ot directly imposed, as in 
other feedback schemes (e .g ISpringel fc Hernquistl 120031 ; 
Oppcnhcimer fc Pavel 120061 ). but occur naturally due to 
thermal pressure exerted on gas in our formulation of su- 
pernova feedback. 
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Figure 7. The angular momentum distributions of stars (blue), 
all gas (red) and HI gas (red dashed) within the virial radius 
at z = and gas which is blown out of the galaxy throughout 
its evolution (black). A bias for low angular momentum to be 
expelled is evident, as is the fact that a large amount of gas is 
expelled. 



The general shape of the outflows histogram is similar to 
the star formation rate, including a particularly prominent 
peak in outflows during the coalescence of the last major 
merger event at time ~ 7 Gyrs (z ~ 0.8). This correspon- 
dence between a merger induced starburst and subsequent 
outflows is discussed in Section [4.41 

The ratio of the outflows from the central galaxy to 
the star formation rate is called the "mass loading" or the 
"loading factor". The log of the ratio of the outflows from 
the central galaxy to the star formation rate is shown in 
Figure U (squares). The mean mass loading from the cen- 
tral galaxy throughout the simulated galaxy's evolution is 
2.3, with a significant increase associated with the merger at 
~ 7 Gyrs. We will see that a high mass loading is a crucial 
ingredient that determines the final angular momentum con- 
tent of our simulated galaxy. In Figure 2] we also show the 
mass loading of the outflows from the major satellite, Sat2 
(crosses), which merges to the central galaxy at z ~ 1. This 
highlights that outflows also occur from the satellites which 
are massive enough to from stars. From here on, we anal- 
yse properties of outflows from the central galaxy, or main 
progenitor, at each time and simply refer to "outflows" . 

In Figure [5] we plot the cumulative mass of such iden- 
tified outflows as a function of circular velocity at the time 
of expulsion. It is apparent that significant amounts of gas 
is expelled while the galaxy is still relatively low mass, but 
that the galaxy continues to accrete material and grow. To 
highlight this, we plot the evolution of the ratio of outflows 
to the total mass of the galaxy (M/ rac ) in Figure [6] This 
plot shows a peak in the fraction of mass being ejected at 
early times, between 2 — 3 Gyrs, when the potential of the 
proto galaxy is relatively low and the outflow rate is high. 
A subsequent steady decline in Mf rac is halted only during 
the last merger at ~ 7 Gyrs, during which large outflows 
occur. 
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Figure 8. The final distribution of j = rVc(r) for "observable" 
baryons in the simulation, which consists of stars and HI gas 
(blue). In red we plot the distribution of all baryons which would 
have ended in the galaxy if outflows were not present, i.e. it is 
the distribution of stars, HI gas and warm/hot gas within the 
virial radius at z = 0, combined with outflows whose angular mo- 
mentum is determined at the timestep prior to being heated and 
expelled. An analytic fit, as determined in Bullock et al. (2001), 
to the typical distribution for dark matter halos is overplotted 
(green). 
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Figure 9. The evolution of the mean angular momentum Q z ) of 
baryons (red) and dark matter (green) which are within the virial 
radius at z = 0, compared with the mean angular momentum of 
gas which is being ejected during the galaxy's evolution (squares). 
Also shown is the mean angular momentum of gas which is be- 
ing accreted at each time (blue line). Material which collapses 
early into the galaxy has relatively low angular momentum, and 
is available for being expelled, and this is replaced by inflowing 
gas which has progressively higher angular momentum. 



4 SELECTIVE REMOVAL OF LOW ANGULAR 
MOMENTUM IN OUTFLOWS 

As a consequence of the strong outflows 70% of the baryons 
are expelled by winds during the galaxy's evolution, and 
there is a significant bias toward low angular momentum in 
the expelled gas. Figure [7| compares the j z distribution of 
gas which has been ejected in outflows (black line) to the 
j z distribution of all baryons within the virial radius of the 
galaxy at z = 0, showing gas (red), HI gas (red dashed), and 
stars (blue). Here, j z is the component of the angular mo- 
mentum whose vector is perpendicular to the disc. The an- 
gular momentum of the outflowing gas is determined at the 
timestep prior to its heating and ejection from the central 
galaxy. The shape of the angular momentum distribution of 
the ejected material differs strongly from the material which 
makes up the galaxy at z = 0. The outflow comprises a large 
low angular momentum peak. The outflows also display a 
significant tail of negative angular momentum material. By 
contrast, the warm/hot gas in the galaxy at z — contains 
high angular momentum material, while the cold (HI) gas 
and stars have very little negative angular momentum ma- 
terial and relatively little low angular momentum material. 
We note that a very similar result occurs when total angular 
momentum is plotted, which we omit for brevity. 

In order to make a more direct link to observations, 
Figure [8] shows the angular momentum distribution, j/jtot, 
for "observable" baryons, defined as a combination of stars 
and c old (HI) gas in order to mimic Ivan den Bosch et al.l 
(2001), who plot such a distribution for several observed 
galaxies. Here, jtot is the mean angular momentum of these 
"observable" baryons, and the angular momentum is derived 
using j = rV c (r), where r is the radius and V c the circular 
velocity. We also plot the distribution of all baryons as they 
would be in the absence of outflows (red line). That is to say, 
the distribution of outflows, as measured just prior to the 



time of ejection, plus all baryons within the virial radius at 
z — 0. This distribution is also normalised by the same jtot 
of the disc baryons, in order to allow direct comparison with 
the blue line which represents the distribution of observable 
baryons. A very different distribution is apparent between 
what is observed at z = 0, and the total baryonic component 
in absence of outflows. In particular a large amount of low 
angular momentum gas is present due to the expelled low 
angular momentum gas, as well as a high angular momentum 
tail of warm/hot gas. 

Also shown on this plot is the distribution of angular 
momentum, normalised by the mean angular momentum, for 
a typi cal dark matte r halo, as characterised in|Bullock et al.l 
(2001) and used in Ivan den Bosch et al.l l|200ll ). The total 
baryons which would have been in the galaxy in the ab- 
sence of outflows (red line) have a very similar shape to this 
distribution, although we note that the normalization, jtot, 
that we used was the mean angular momentum of the disc 
baryons. Normalising by the mean of the total angular mo- 
mentum of baryons which would have been in the galaxy 
at z = would shift the distribution to the left. Neverthe- 
less, the shape of the distribution of all baryons which would 
have been in the galaxy in the absence of outflows, and that 
of dark matter halos, are remarkably similar, highlighting 
the role that outflows play in causing the distribution of 
observed baryons to differ so sharply from their host halos. 

We now identify the processes which result in the strong 
bias toward low angular momentum of outflows that we have 
identified. In combination they result in the ejection of low 
angular momentum material and a significant alteration of 
the angular distribution of baryons, and in particular the 
distribution of stars and cold (HI) gas within galaxies com- 
pared to their host dark matter halos. 
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4.1 Early accreted material has low angular 
momentum. 

Material which falls into the central regions of proto-galaxies 
at the earliest times has been subjected to the least amount 
of torquing from the large scale structure (|Rvdenl 1 19881 ; 
iQuinn fc Binn^vll 1992ft . Material falling in at later times has 
higher angular momentum. This results in early outflows 
having less angular momentum than later outflows, and in- 
troduces a bias whereby low angular momentum material 
is preferentially expelled. Figure [9] shows the evolution of 
j z of outflows (squares), calculated at the time of ejection, 
as well as the inflows (blue line) at a corresponding time. 
Also shown is the evolution of j z for the baryons (red line) 
and dark matter (green line) which end up in the galaxy at 
z — 0. A similar plot results when total angular momen- 
tum is plotted in this manner. The central galaxy, or main 
progenitor, is used as a reference frame at all times, and as 
we stated above, outflows are those ejected from the central 
galaxy, while inflows, described below, are also to the central 
galaxy. 

After spinning up, the dark matter material which ends 
in the galaxy at z — essentially retains its angular mo- 
mentum. The baryons which end up in the galaxy at z = 
have also retained their angular momentum. At the early 
times, it is low angular momentum material that has been 
accreted and is hence available for blow out and this is re- 
flected in the very low angular momentum of material which 
is ejected at early times. Only gas which has been accreted 
can be ejected, so a large bias is clear between ejected gas 
and the gas which will end up in the galaxy at z = 0. The 
low star formation rates and relatively high blowout rates 
(i.e. the high mass loading, see Figure[4]| means that a large 
fraction of low angular momentum material is not forming 
stars. As outer, high angular momentum shells of gas are 
accreted, they make up a larger fraction of the angular mo- 
mentum distribution of gas than they would if none of the 
earlier accreted gas were expelled. This naturally alters the 
distribution of gas angular momentum compared with dark 
matter. Low star formation efficiency and high gas fractions 
mean that blown out gas has preferentially low angular mo- 
mentum. 

The blue line traces the angular momentum of inflowing 
gas at each time. We define inflowing material as gas which 
has reached 57?* for the first time, in order to make a fair 
comparison with the outflowing material as we have defined 
it (recall that 7?, is the radial distance to the furthest young 
star from the centre of the galaxy at each time) . Clearly, in- 
flowing gas has significantly higher angular momentum than 
the outflows at all times. At late times, the inflowing mate- 
rial has higher angular momentum than the mean of all gas 
which is within the virial radius at z = 0, reflecting the fact 
that later accreted material has high angular momentum. 

Lower potential wells at early times favour higher 
blowout rates, at a time when lower angular momentum ma- 
terial is being accreted. 

4.2 Extended reservoir of high angular 
momentum material. 

Star formation only occurs in the inner regions of the 
galaxy where the gas reaches sufficiently high densities. Ex- 
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Figure 10. The face-on B-band surface brightness map of the 
galaxy at Z ~ 1.2, over plotted by an HI contour map. The HI con- 
tours range from 10 19 cm — 2 to 10 22 em -2 in steps of 10°' 5 cm -2 . 
The extent of the star forming region is reflected in the B-band 
map, while the extended nature of the HI is evident. 

tended HI gas discs are found arou nd local isolated low 
mass galaxies (|Broeils fc Rhed 1 19971 ). Such low mass, iso- 
lated galaxies have gas dominating their baryonic mass frac- 
tions, with values typically above 70%, a nd as high as 95% 
l|Schombert et al.ll200ll ; iGeha et ajj|2006ft Ther e is evidence 
that extended HI discs occur in i ntermediate llPuech et all 
l2010ft and high redshift galaxies (<Daddi et al.ll2010ft . Fur^ 
ther, there is evidence for the existence o f warm /hot gas in 
the halos of disc galaxies (|Spitzerl 1 19561 ; ISpitzer fc Jenkinsl 
Il975l ; ISembach et al.ll2003ft . In our simulations, gas beyond 
the star formation region serves as a repository of high an- 
gular momentum material. Typically, between 30% and 40% 
of gas within the virial radius of our simulated galaxy lies 
within the star forming region at a given time. Figure [10] 
shows the face-on B-band surface brightness map of the 
galaxy at z ~ 1.2, over plotted by an HI contour map. 
The extent of the star forming region, from which outflows 
are driven, is reflected in the B-band map, while the ex- 
tended nature of the HI is evident. At this time, 47% of all 
gas is beyond the star forming region. Gas within the star 
forming region has average specific angular momentum of 
jz m ean = 66 kpc kms -1 while gas beyond the star forming 
region has jzmean = 131 kpc kms 1 . Looking at just the HI 
gas at this time, 20% of HI is beyond the star forming region, 
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Figure 11. Left Panel: The background is a density map of all gas in the simulated galaxy at z = 0.5, where the disc is oriented edge 
on. We have taken a 200x200 kpc box. The dashed line marks the virial radius. The arrows are velocity vectors, indicating the direction 
of outflows, with the size of the arrows related magnitudes of velocities as indicated above the panel. Middle Panel: The gas which is 
feeding star formation, i.e. we include only gas at z = 0.5 which will subsequently form stars by z = 0. We use the same scale and again 
indicate the virial radius. The gas feeding the star formation is accreted largely from low angles from the plane of the disc. Right Panel: 
Here we plot outflowing gas, again with velocity vectors demonstrating that outflows are strongly directed perpendicular to the disc. 



and has jzmean = 145 kpc kms -1 . Gas which gains energy 
directly from supernovae thus has relatively low angular mo- 
mentum, so simply modelling this gas which is within the 
star forming regions as being blown out will preferentially 
remove low angular momentum gas. 



4.3 Outflows perpendicular to the disc. 

Our simulations, as well as observations of outflows, indi- 
cate that as gas is expelled, it can also entrain gas which 
is in outer r egions of the galaxy, and sweep it out with 



is m outer r egions ot th e galaxy, and sweep it out with 
the outflows ([Stewart et al.ll2000l ; ISchwartz fc Martini liool 



IVeilleux et all T2005). The relative importance of these two 
modes of outflow are not well constrained. Direct expul- 
sion from inner regions appe ars to be the dominant mode 
l|Bland- Hawthorn et al J 12007) . and is also the main mech- 
anism of gas removal in our simulations. Here we show 
that t hat, consistent with observations ([ Heckma n et al.l 
1987 1; iBland fc Tullvl ll98Sl; IShopbell fc Bland-Hawthorn! 



199a : IVeilleux fc Rupkell2002l ; [Bland-Hawthorn et al 



and other theoretical models ( Mac Low fc Ferraral 



2007) 



1999 



iPieri et al.l 120071 ). the outflows from the central regions are 
preferentially perpendicular to the plane of the disc. 

Figure [TT] shows the direction of gas flows from our 
galaxy at z ~ 0.5, when a minor starburst occurs. The color 
gas density map shows the edge on view of the galaxy, with 
velocity vectors overplotted showing the direction and mag- 
nitude of the gas velocity. The virial radius is marked by 
a dashed white line. The right hand panel shows the out- 
flows, and reveals that they follow the path of least resis- 
tance perpendicular to the disc. These outflows thus prefer- 
entially entrain material in regions which are perpendicular 
to the disc. The highest angular momentum material, which 
is in the extended gas disc surrounding the star forming 
regions, is the least affected by the outflows. We measure 
a low value for the mean angular momentum in the plane 



of the disc of this material which is identified as outflows, 
with jz m ean{put flows) = 20kpc kms -1 , while total angu- 
lar momentum is high, jmean(out flows) = 193kpc kms , 
as the material is getting driven to large radii, perpendic- 
ular to the plane, at high velocities. Here, we have used 



Jmean 



V. 



By contrast, gas which feeds star formation falls on 
the galaxy from a direction which is in the plane of the 
disc, and away from the outflows. The central panel of Fig- 
ure [TT] shows the edge on galaxy, again at z = 0.5, but this 
time only gas which will form stars by z — is plotted, 
again with vectors indicating the velocity of the gas. The 
gas which feeds star formation is clearly shown to flow pri- 
marily from the plane of the disc. This gas has a relatively 
high planar angular momentum, jz m ean(starfeed) = 79kpc 
kms -1 , which is a large fraction of this material's total an- 
gular momentum j rnean (starfeed) — 81kpc kms -1 . This 
contrasts to the outflowing material, but also to gas which 
is in the galaxy at this time and which remains in the 
galaxy as gas at z = 0, i.e. material which is neither be- 
ing ejected nor feeding star formation. Such "retained" gas 
has jzmean (retain) = 67kpc kms -1 , with total angular mo- 
mentum, jmean(retain) = 127kpc kms -1 . 



4.4 What about mergers? 

One of the characteristic features of CDM cosmologies is the 
merger of galaxies. This process is ubiquitous, and essen - 
tially scale free (e.g. ICole et aL l l2000l : lFakhouri et alJlioich . 
meaning that low mass field galaxies which preferentially be- 
come discs dominated will have in general similar merging 
histories as high mass galaxies, which have bulges and may 
be dominated by spheroidal star systems. Mergers are widely 
expected to result in bulges. In earlier simulations of galaxy 
formation which suffered from overcooling, proto-galaxies 
rapidly formed central baryonic concentrations through dy- 
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namical friction (|Navarro fc Whitdl 19941 ) This problem was 
exacerbated in lo w resolution simula tions, which have denser 
central regions (|Maver et alJ 120081 ). As explained above, 
feedback from supernovae prevents a dense bulge forming in 
the central regions of the proto-galaxi es in these current sim - 
ulations. Further, as demonstrated in IStewart et ail (j2009h . 
direct accretion of stars is insignificant in low mass galaxies, 
whose mergers are dominated by gas. The stars which are 
accreted in the mergers, tend to get thrown into high en- 
ergy orbits during the interaction, meaning that the result- 
ing spheroid is not centrally concentrated. A diffuse stellar 
halo rather than centrally concentrated bulge results from the 
scrambled trajectories of stars during the late major merger 
in our simulation. Indeed, at z=0, our simulated galaxy has 
an old (mean age 10.3 Gyrs), diffuse stellar spheroidal com- 
ponent which comprises only ~ 6% of the stellar mass of the 
simulated galaxy, even though it has no bulge. This is con- 
sistent with observed low mass galaxies in all environments 
which have old, dim spheroidal st ellar populations in ad- 
dition to their main body of stars llVansevicius et al.l |2004| ; 
iHidalgo et al. l l2009l;IStinson et alj|2009h. 

Simulations, (Governa to et alJ |2009; Hopk ins et al.l 

2009) showed that, when mergers are gas rich and feedback 
is present, bulges can be as little as 10% of the resulting 
galaxy. In such gas rich mergers, no strong bar forms to 
drive gas to the central regions. Our simulated galaxy has 
a rich merger history including a merger with mass ratio of 
3:1 which has a final coalescence z ~ 0.8. All mergers are gas 
rich, with gas to stellar mass ratio above one for all progen- 
itors of our simulated galaxy. Several reasons ensure that a 
bulge does not form during these mergers. Firstly, only the 
inner cold gas is subject to being driven to the centre dur- 
ing the mergers. The surrounding gas is not affected, but we 
know that inner gas has preferentiall y low angular momen - 
tum. Secondly, in our simulation as in iHopkins etal](|2009l ). 
a bar does not form in gas dominated mergers, limiting the 
amount of gas driven toward the centre. 

Thirdly, and most importantly, the feedback from the 
small fraction of such gas which forms stars in a burst expels 
a large fraction of the gas which has lost its angular momen- 
tum and has been driven toward the centre of the galaxy, 
as reflected in the high mass loading during this merging 
epoch, shown in Figure 2] This is clear from the dip in the 
line traced by squares in Figure [9l which indicates that the 
material which loses angular momentum material is removed 
during the merger. Star formation is primarily occurring in 
the very inner regions of the merging galaxies, and hence it 
is coupled to the very material that has lost angular momen- 
tum which is being blown out. The effective coupling of SF 
with large scale supernova driven outflows during m erger in- 
duced star bursts is well supporte d observationally l|Martinl 
ll999l ; lBland-Hawthorn et alj|2007ft and is the extra ingredi- 
ent in our models which allow us to form bulgeless merger 
remnants. 



creates pressure, and is modelled to undergo a blastwave 
phase which drives galactic winds. We emphasise that nei- 
ther degree nor direction of galactic winds are enforced a 
priori within our supernova feedback scheme. It is necessary 
for the spatial resolution to be significantly le ss than the 
scale length of the disc to capture these effects (|Coh'n et al.l 
2010). More gas is expelled from our simulated galaxies than 
forms stars by z = 0, up to 70% of the original cosmic frac- 
tion. Yet, similar to observed isolated dwarfs dGeha et al.l 
2006), they remain gas dominated at all times, due to the 
continued cooling of gas into the disc and the low star forma- 
tion rate. We have previously shown that our model results 
in simulated dwarf galaxies which match multiple properties 
of observed dwarf galaxie s (|Governato et al ] |20ld ; IOh" et al.l 
l2010l; iBrooks et al.ll2O10h . Future and ongoing work will in- 
clude metal lines cooling at all temperatures, H2 cooling 
and its form ation and destruct i on, and SF linked to the local 
H2 fraction jBooth et aL I l2007l ; iRobertson fc Kravtsovll2008l ; 
iGnedin et afl l2009V Preliminary results show that the high 
density SF threshold adopted here closely mimics the dis- 
tribution of SF in simulations where SF is allowed only in 
gas with a high fraction of H2 (Christensen et al, in prep) . 
The additional cooling by metals and H2 is counterbalanced 
by an increase in the SN efficiency (Pontzen et al, in prep). 
None of these details of localised star formation affect the 
galaxy scale processes which are highlighted in this current 
study. Here, we have outlined the reasons that the galactic 
winds in our simulated galaxies are strongly biased toward 
expelling low angular momentum gas. We identified the in- 
terlinked processes responsible for such bias, which naturally 
occur when feedback is effectively modelled in a hierarchi- 
cally assembled galaxy in a CDM Universe, namely: 

• Low angular momentum material is accreted and 
rapidly expelled early, when SF peaks 

• A low potential at early times favours early gas ejection. 

• Later accreted material has higher angular momentum 
and forms a disc. 



These ideas are similarly outlined bv lDekel fc Silkl l|l98rj ); 
iBinnev et al.l (|200ll ) , where a notion of expelling gas early is 
invoked to explain the angular momentum of disc galaxies. 
Early ejec tion of material is also postu lated in recent 
models of iDutton fc van den Boschl l|2009l ) to explain the 
mass distribution of disc galaxies. 

• A reservoir of high angular momentum material exists 
beyond the star forming regions. 

Di sc galaxies have extended HI discs jBroeils fc Rheel 



1997h. as well as warm / hot gas in the halo llSpitzerlll956r 
Spitzer fc Jenkins! 11973 : ISembach et~aH l2003h . These act 



as reservoirs of high angular momentum gas, which are 
beyond the star formation radius and thus not subject to 
being directly expelled from the galaxy. 



5 DISCUSSION 

In our simulated dwarf galaxies, we are able to resolve the 
process of energy injected from multiple overlapping super- 
novae, which occur in star forming regions coupled to dense 
gas. The thermal energy injected into the surrounding gas 



• Outflows occur perpendicular to the disc, entraining 
relatively low angular momentum material and preventing 
infall of material from regions above and below the disc 
plane. 

Observations strongly indicate that galactic winds expel 
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mate r ial perpendicular to galactic discs dHeckman et "al 



19871 ; i Shopbell fc Bland-Hawthornlll998l ; IVeilleux fc Rupkc. 
20021; bland-Hawthorn et al.|l2007h. as predic ted by theory 
|Mac Low fc Ferraral 1 19991 ; IPieri et all 120071 '). Gas which 
builds the disc at late times is accreted from regions that 
are relatively aligned to the disc plane. 

• Mergers or instabilities which cause gas to lose angular 
momentum, trigger star bursts which expel most of this low 
angular momentum gas. 

Star burst galaxies are generally associated with 
mergers, and are ob s erved to have significant outflows 



dHeckman et al.l Il990l: | Griffiths et al. 1 12000| ; IVeilleux et al 



2005) in local galax i es an d at high redshift (jPettini et al 
20001 ; iNestor et al] l20ld) , pre dominately driven from 
the inner starburst regio ns dChevalier fc Clegd 1 19851 ; 
IStrickland fc Heckmar] 120071 ). ensuring that it is the very 
material that has lost angular momentum during the 
merger which is being expelled. High velocity outflows 
from pos t-starburst galaxies are also observed at high 
redshift (ITremonti et al.l 120071 ). The high mass loading 
in these events prevents large amounts of stars being 
formed from this material which had been driven toward 
the galaxy's inner regions. Feedback from central black 
holes/ Active Galactic Nuclei (AGN) would also aid in 
expelling material which has lost angular momentum and 
has been driven to the central regions. There is evidence to 
sugge st that such feedback is directed perp endicular to the 
disc l|lrwin fc Sofuel 1 19921 ; ISu et all |2010T ). enhancing the 
processes outlined in this study. 

• A diffuse stellar halo rather than a centrally concen- 
trated bulge results from the scrambled trajectories of stars 
during the late major merger in our simulation. 

The stellar rather than gas component of mergers between 
progenitors which are massive enough to contain stars them- 
selves, but which do not themselves have bulges, result in a 
diffuse halo component. 

By expelling large amounts of preferentially low angu- 
lar momentum gas from the galaxy, particularly at early 
times and during merger events, the shape of the angular 
momentum distribution of baryons is altered significantly 
from that of their parent dark matter halos. This allows the 
formation of bulgeless disc galaxies. Future work is needed to 
see if the same processes (possibly in conjunction with black 
hole/ AGN related feedback which also emanates from cen- 
tral regions and thus expels low angular momentum gas) is 
able to explain the suppression of bulges in more massive ha- 
los and the f ormation of bright disc galaxies with small clas - 
sical bulges l|Kormendv et al.ll2010l ; [Peebles fc Nusserll2010l ). 
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